Effects of finite probe size on self-affine roughness measurements 
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The roughness of fracture surfaces has been shown to exhibit self-afhne scale invariance for a wide 
variety of materials and loading conditions. The range of scales over which this regime extends 
remains a matter of debate, together with the universality of the associated exponents. The topog- 
raphy of these surfaces is however often investigated with a contact probe that is larger than the 
micro-structure. In this case, we show that the correlation function of the roughness and the corre- 
sponding Hurst exponent can only be measured down to a length scale Axc which depends on the 
probe size R, on and on the surface topothesy I, and exhibit spurious behavior at smaller scales. 
First, we derive the dependence of Aa;;. on these parameters from a simple scaling argument. Then 
we study this dependence numerically and verify our theoretical prediction. Finally, we establish 
the relevance of this analysis from AFM measurements on an experimental glass fracture surface 
and provide a metrological procedure for roughness measurements. 
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Fractography has long been a fruitful method to deci- 
pher the microscopic mechanisms of fracture in hetero- 
geneous materials (see e.g.,|Tll2] for reviews). Since an 
early paper by Mandelbrot et al. a growing interest 
has developed in the understanding of the fractality of 
crack surfaces. One of the most prominent features of 
these surfaces is the self-afRne structure of their rough- 
ness in various materials, ranging from ductile metallic 
alloys to brittle mortars or glasses (4l l5l [6l Ffl [8l [9l ITOl ITT] . 
A proper understanding and experimental characteriza- 
tion of this feature has both applied and fundamental 
stakes: on the one hand, it has recently led to a method 
to extract the direction of crack propagation from the 
analysis of roughness, which is of high practical impor- 
tance [12]. On the other hand, the possible universality 
underlying this critical-like behavior remains to be as- 
sessed [13] , together with the physical origin of the intrin- 
sic cutoffs in the self-affine regimes [5J [T2] . These scaling 
properties have motivated the development of numerous 
models O US US HH [H US HU] ■ 

Recent advances in scanning probe microscopy (SPM), 
especially the development of atomic force microscopy 
(AFM), have provided large amounts of experimental 
data drawing a rather complicated picture of the situ- 
ation: the reported exponents depend on the direction of 
analysis, the material and the resolution. Furthermore, 
the characteristic length scale of the self-affine regimes, 
or topothesy, is rarely reported or unphysical. This dis- 
crepancy in the reported exponents can originate from 
the existence of several universality classes and actual 
differences in the underlying physical mecanisms. How- 
ever, it might also partly result from experimental biases. 
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The list of possible biases is long and instrument depen- 
dent, i.e. environment noise, feedback loops, biases in the 
estimators [TB] , etc. Here, we focus on the subtle bias in- 
duced by the non-linear smoothing effect of the finite-size 
contact probe used in topography measurements of self- 
affine surfaces [2Tj. 

In order to quantitatively estimate this bias, we in- 
vestigate the consequences of using model probes of var- 
ious sizes and shapes on ideal self-affine surfaces. As 
SPM topography measurements usually proceed line by 
line, we will study the one-dimensional problem of prob- 
ing a self-affine profile h (x) with a rigid tip defined by 

/]^' (x) = 2i||— T, where R is the characteristic size of the 
apex and 7 a shape parameter (7 > 2). This choice of 
shape is motivated by the possible flattening of real tips 
during scans. The roughness features of the profile h 
are extracted from the height-height correlation function 
defined as 



Ah (Ax) = d {{h {x + Aa;) - h {x)Y), 



(1) 



The profile h{x) is self-affine if it is invariant under the 
transformation 
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where C is the roughness or Hurst exponent. It follows 
that there is a physical length scale I, called the topothesy, 
for which the local slope of the profile becomes one. It 
verifies 

C 

i.e Ah = Ax^l^-'^ (3) 
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Physically, this parameter provides information about 
the amplitude of the roughness. In particular, it follows 
from Eq. ^ that vertically dilating the profile affects the 
topothesy as 



lih I — > fj, i-c I 



(4) 
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Measuring such a profile with a finite-size SPM tip 
amounts to finding the height, h{x) of the apex as the 
tip first touches the surface when approached from above, 
for each value of depicted in Fig. [T] This occurs at 
some position xq such that 

h{x) = Sji [h] (x) = max^' {h (x') - gnix- x')) (5) 

where gn (x) describes the shape of the tip, R its size, 
and Sfj is a smoothing operator acting on profiles. Sn 
is non- linear (Sji [fih] ^ (J-Sr [h] ) and implies partial loss 
of information. Interestingly, for tip shapes gfj — /jj'' , 
we can still derive a scaling relation for Sr [hi. Writing 
Eq. (|5| for a profile scaled according to Eq. Q yields 



S]i [fih] = maxj,/ I ^h{x') — 



Sr [fJ.h] = M < maxj;/ h{x') 



\x - x'p 
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2 l^fi'^R] 

SR[nh] = fiSR'[h] with R' = fi^R (6) 
This relationship states that the bias induced by scan- 
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The root mean square of this quantity, Ah, scales like 
Ax, with a prefactor that we conjecture to be propor- 

tional to ^g- (Axc)- On the other hand, the very large 
scale structure of the profile is not affected, and we ex- 
pect Ah ~ Ax'^ for Aa; ^ Axc- We can now evaluate 
how Axc scales with R, ^ and I by assuming that at Axc, 
the height of the tip equates the height-height correlation: 



/jj^l iAxc)-Ah{Ax,) 



(7) 



Using Eq. ^ and the expression for the chosen tip shape 
yields 



2R'r-^ 



(8) 



Solving for Axc and dropping the numerical prefactor 
gives the approximate scaling relation 



Axc {R'^-H^-'^)-'-' 



(9) 



In order to assess the validity of this relation, we 
have performed numerical simulations of self-affine pro- 
files with exponents C ranging from 0.2 to 0.9 and scaled 
the topothesy I in the range 10~'^/o to lO^Zg according 
to Eq. Q. Due to its high representativeness, we focus 
on the parabolic case 7 = 2 with tip radius of curvature 
R ranging from 10~^lo to 10^ Iq, where Iq is the spatial 
sampling interval. In the following, we will drop Iq in our 
notations, and consider all length scales in units ^o- The 
self-affine profiles have been obtained using the spectral 
method described in J22} ; they contain 20000 data points 
each. For all values of these parameters, we compute 
the height-height correlation function for the original and 
the smoothed profiles. An example is shown in the main 
panel of Fig. ([2|. The self-affine exponent and topothesy 
of the simulated profiles are evaluated according to the 
following method. We define the function A by 



FIG. 1: Schematic of the small scale behavior of the smoothed 
profile. For two closeby points, the probe (black) touches a 
unique distant peak on the profile (blue). We have drawn the 
corresponding portion of the measured profile h (dotted red) 



ning a profile with a given tip-radius is equivalent to 
that obtained by scanning the vertically expanded pro- 
file with a thinner tip. Now we turn to the investi- 
gation of the influence of the local probe size on the 
height-height correlation function. Qualitatively, apply- 
ing S dampens the roughness details at length scales 
smaller than a characteristic length Axc- Let us consider 
two nearby points xi and X2 close enough to a point 

Xq such that h{xi^2) — h{xo) — (a^i,2 — a;o), as de- 
picted in Fig. [T] Setting X2 = Xi + Ax, for values of 
Ax much smaller than Axc, we get h(x2) — h(xi) 

Ax) ~ (xi - xo) 



xo 



%(xi-xo)Ax. 
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(10) 



It follows from Eq. ^ that A (Ax, C) = I for all Ax. 
We evaluate A (AX, rj) for a uniform random variable 
AX. C, is then defined as the value for which the vari- 
ance of A reaches a minimum in rj, and the average value 
of A (AX, C) provides an estimator of I based on its fun- 
damental invariance property Aft, (/) — I. 

We observe a clear crossover between the small and 
large scale behavior of Ah. At large scales. Ah asymptot- 
ically behaves like Ah. It has to be noted that it reaches 
this regime very slowly, and even though the exponent 
reaches C from above. Ah remains underneath the origi- 
nal correlation function, resulting in underestimation of 
the surface roughness amplitude and overestimation of 
the exponent. At scales smaller than the crossover length 
scale Axc, Ah bends down and reaches the expected lin- 
ear scaling after a long transient. In order to extract the 
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position of the cutoff Axc, we extract the local behavior 
of A/i around Ax = Iq by fitting a line 



logio Ah = a{R, C, logio Ax + (3{R, C, I) 



(11) 



through the first two points of A/i: a is essentially 
the log-derivative of h at the resolution Iq (i.e. at the 
smallest accessible value of Ax). We have represented 
in Fig. [sj^Top) the behavior of a and /3 as a function 
of the tip radius R for C = 0.33 and / = 100. Two 
regimes are clearly observed. For values of R smaller 
than a cutoff -Ro(Cj 0; the height-height correlation func- 
tion is not affected by the smoothing effect. From 
Eq. ([3]), we get that in this regime, a{RXJ) = C ^^nd 
/3 (i?, C, = (1 - C) logio I when R < i?o(C, 0- For values 
of R much larger than Ro, a converges towards 1, as ex- 
pected, whereas P essentially behaves like a power law in 
R. In order to extract the scaling of the threshold Rq, we 
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FIG. 2: Main panel: Height-height correlation function eval- 
uated for i) a self-affine profile with = 0.33 and I = 100 
(•), ii) the same profile smoothed with a tip of size R = 7.8 
(■). Inset: Portion of the corresponding original (blue) and 
smoothed (red) profiles. 



It is interesting to notice that since C — 1 is smaller than 
1, the effect of the tip appears for smaller radii when 
the topothesy is larger. This scaling is a consequence of 
the equivalence expressed in Eq. (|6|. The existence of 
this cut-off is an effect of the resolution Iq of our mea- 
surement. As a matter of fact, restoring the sampling 
length scale yields Ax^ (Ro) = lo- This means that for 
radii smaller than Rq, the cut-off length Axc becomes 
smaller than the resolution and hence disappears from 
the resolved portion of the correlation function. In the 
asymptotic regime where i? 3> i?o, the small scale expo- 
nent a goes to 1, as can be seen in Fig.jsjTop). Thus the 
asymptotic form of the correlation function for R^ Rq, 
or equivalently Axc ^ Ax, reads: 



logio Ah = logio Ax 



l-C, I 
^ logio ^ 



Solving Ah{Axc) = Ah{Axc) yields Eq. ( 
above for the expression of Axc. This resu 
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fit /3 = a {(, I) logio R + b{C,l)- The function a {(, I), ap- 
pears to depend only marginally on I, which we neglect. 
The dependence of a on C is depicted in Fig. |3] (Bottom 
left). It can be accurately modeled by a = — We 
now turn to the scaling of b, which is well described by 
b = u logio ^ + No significant dependence on C, is found 
for V, and this number is of the order of —0.1, which 
is negligible and will be dropped henceforth. We find 
u = ^E^, as shown in Fig. 3 (Bottom right). Altogether, 
the behavior of /? conveniently writes 

/3 = (^) logic ^ when R » Rq (12) 

We can now extract Rq by intersecting the small and 
large scale behaviors of /3, which reads, in units Iq 

(1 -0 logio logio (13) 

This yields 

logio i?o = (C-l) logio ^ (14) 



FIG. 3: Top: Apparent parameters a and /3 as a function of 
logj^Q R for ( — 0.33 and / = 100. Bottom left: Parameter 
a as a function of exponent C,. The line represents —jz^- 
Bottom right: Parameter m as a function of exponent 
The line represents ^E^- 

our asymptotic conjecture in the parabolic case since 

r//"'^' At 

logio ^ (Ax,) = logio ^=/3 (16) 

To check the relevance of this analysis, we have per- 
formed AFM (Veeco Nanoscope Dimension V) measure- 
ments on fused silica glass (Corning 7980) fracture sur- 
faces. We break a DCDC sample in the stress-corrosion 
regime, in conditions detailed in [8l. Then we scan the 
topography of the fracture surface with the AFM in tap- 
ping mode (2 x 2^m^,512 x 512 pixels). The height- 
height correlation function of the data parrallel to the 
crack front is shown in Fig. |4] Two scaling regimes are 
observed: a 0.8 exponent at small scales, and a 0.2 ex- 
ponent at large scales. Once this is done, we acquire an 
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image of the tip of the AFM with a silicon tip character- 
izing grating (TGTOl, NT-MDT, Russia). By fitting /jj'' 
on this measurement, we evaluate Rxip = 40 ± 5nm and 
a shape parameter 7 = 4± 0.5, which is typical for AFM 
tips on hard substrates. Then, through trial and error. 
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FIG. 4: Experimental height-height correlation functions ob- 
tained for a silica glass fracture surface scanned with an AFM 
with tip radius R — 40nm. Original and smoothed height- 
height correlation functions obtained for a numerical profile 
with = 0.13 and a 40nm tip for different values of 7. 

we select a numerical profile with ^ = 0.13, which is close 
to the classical logarithmic behavior [T^. We smooth it 
using the measured value of Rxip, and different values 
of 7. The resulting correlation functions for this value 
of Ci superimposed to the experimental ones in Fig. |4] 



reproduce quantitatively our measurements for 7 = 3.5, 
which is compatible with the measured shape of the tip. 

Altogether, these results strongly suggest that the 
small scale self-affine regime with a large exponent and 
the associated cut-off we observe on silica glass fracture 
surfaces are due to the tip-smoothing effect. We conjec- 
ture that earlier measurements using contact probes on 
different materials might exhibit such bias. The numer- 
ical procedure we proposed provides a sharp metrologi- 
cal constraint on roughness characterization and allows 
for accuracy checks on earlier results. Further work is 
needed to understand the slow convergence of the corre- 
lation function at large scales and reconstruct the statis- 
tical properties of the original roughness from a smoothed 
profile. 

To conclude, we have shown that the use of a finite-size 
contact probe to scan a self-affine surface significantly af- 
fects its measured small scale roughness properties. In 
particular, we have derived theoretically and checked nu- 
merically the dependence of the length scale at which this 
effect starts on the surface parameters and the asymp- 
totic small scale behavior of the height-height correla- 
tion function. Finally, we have performed experiments 
confirming the relevance of this analysis regarding AFM 
metrology and provided insight into the origin of previ- 
ously reported large Hurst exponents at small scales. 
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